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Abstract 

In this work stability of polygonal configurations on a plane and sphere 
is investigated. The conditions of linear stability are obtained. A nonlinear 
analysis of the problem is made with the help of Birkhoff normalization. Some 
problems are also formulated. 

1 A linear stability of Thomson's configurations 
of vortices on a sphere 

Motion equations of a system of iV point vortices on a sphere can be presented in 
the Hamiltonian form pQ with the Poisson bracket 

{<Pu cos0k} = (1) 

and the Hamiltonian 

N 



i/ = -i f yj'r i r t i„(4tfsm 2 f) 



i,fc=l 



(2) 



cos 7^ = cos 9i cos 9k + sin 9i sin 9k cos(c^ — tpi 
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where fi and Tj are spherical coordinates and the intensity of the z-th vortex, 
7,fc is an angle between the z-th and the k-th vortices, R is a radius of the sphere the 
prime denotes that i ^ k. Here and further all indices values from 1 to N. Motion 
equations of such a system have the form 

N 

■ 1 sin Qi sm(f k - if A 

U k = —. — Fv7 / , 1 i 



AuR 2 ^—^ 1 1 — cos 7^ 

(3) 



N 



sm 9 k f k 



1 ^-v -n s i n @k cos Qi — cos 9 k sin 9i cos(fk — fi) 



AnR 2 1 - cos 7^ 

In the case of equal intensities Tj the system is invariant with respect to the 
discrete group of permutations of vortices, and therefore allows different symmetric 
partial solutions. Let us consider partial solutions, which are analogous to the well 
known Thomson's configurations on a plane. Then the vortices are located on the 
same latitude 9q in the vertices of a regular iV-gon and rotate around its center with 
the angular velocity uj: 

r, /„s 27T,. „ s .„ . T(N — 1) COt 0n , 

< = O = const, fi = u{e )t + —{t-l), u{9 ) = \ m ■ a ■ 4 

JM 4nR z sm t/ 

Let us consider stability of these partial solutions in linear approximation. The 
study of stability of Thomson's configurations of vortices on a plane in linear ap- 
proximation was carried out by J. Thomson [2J. He has shown that they are stable in 
linear approximation for N < 7, and unstable for N > 7. Let us take variations 59i 
and Sfi as 

9 l = 9 + 59 l , fi = u(9 )t + — (i - 1) + Sfi . (5) 
The linearized equations for 59i and Sfi are autonomous and have the form: 



T 1 

4ttR z 2 sm 6*0 



r 1 A c „ T(N -1)1 + cos 2 9 
where A is an (N x iV)-matrix with the elements 



(6) 



sm9 6f k = — — j o , a AkiS'&i — 7 - T - — 59 k 

AnR 2 2 sin 2 9 AnR 2 sm 2 9 



N-l 1 1 
A ki = Ski^2—^ (1 ~ hi) . 2 nf U T ' 

^— ' sm Tfin sm 4?(k — i) 

By differentiating with respect to time, the equations © can be reduced to a system 
of N second-order equations in the following form: 

Sf k = I 8lvR2 T sin 2 g Q J (4y ~ 2(iV - 1)(1 + cos 2 9 )S kj ) AjiSfi (8) 
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(similar equations are obtained for 9k under elimination of tpk). 
A solution of the equation (jSJ) has the form 



(m)n m t 



where the constants C^, Q m are expressed through the eigenvalues and eigenvec- 
tors of the matrix A^. As in the case of a plane the elements of Aki depend 
only on differences (k — i), thus the matrix A can be diagonalized by the Fourier 
transformation, and the eigenvalues of the matrix can be presented as 

r Ah 

A m = 2m{N - m) , m = 0, 1, . . . , [— J . (9) 

Each eigenvalue corresponds to two eigenvectors 4>f^ = cos f^fp*) and = 

sin (^p«), i = 1, • • • , N. To the values m = and m = y (in case of even AO, 
there corresponds a unique eigenvector cf> \ jH] . 
Using (jHJ), ©, we find frequencies fl m 



^ m = Q D2 m 2fl V / A m -2(iV-l)(l + cos 2 go 
87rit^ sin Vn 



— m) 



(10) 



y/m(N - m) - (N - 1)(1 + cos 2 6» 



47ri2 2 sin 2 

The presence of the zero frequency (II (J |) at m = corresponds to instability of 
Thomson's configurations in the absolute coordinate system. Indeed, under pertur- 
bations 89i, 5 <pi corresponding to m = the Thomson's configuration transforms 
into Thomson's configuration, which is close to the latter and moves away from 
it linearly with time. Stability of relative motions of the vortices, from which we 
exclude the rotations of the system as solid-state configuration, are determined by 
the remained frequencies. It is necessary for the relative stability the remained fre- 
quencies Q m purely imaginary. It is obvious from (JTUJl that the configurations (@J) 
gradually lose their stability approaching the equator (i. e. with magnification of 9q). 

Table 1 shows the boundaries of varying of the radicand in (fTUj) for various Af 
with the maximum value m(N — m) and also values of 9$, at which the radicand 
changes the sign. 

For N > 7 there is always at least one frequency with a positive real part, hence, 
such configurations are unstable already in the linear approximation. 

For A^ = 6, 5, 4, 3 and 2 and for # G (0, 9q) all frequencies Q m are purely 
imaginary. Thus, under the above-mentioned conditions the specified configurations 
are stable in linear approximation. For #o G [9q, |] these configurations are either 
unstable in linear approximation or consideration of higher order approximations 
(at 9 = 9q) is required for determination of stability. 
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Table 1. 
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2 A nonlinear stability of Thomson's configura- 
tions on a sphere and a plane for AT = 2, . . . , 6 

The work by L. G. Khazin [3] introduced a proof of the Lyapunov stability of Thom- 
son's configurations in nonlinear statement on a plane. However, though the general 
theorem from the stability theory, proved by L. G. Khazin in this work is correct, 
its application to Thomson's configurations on a plane is not appropriate. Zero fre- 
quencies considered by Khazin correspond to the motion of vortices as a whole and 
are not connected with the relative stability of Thomson's configurations. Thus, 
the problem of stability of Thomson's configurations on a plane in strict nonlinear 
statement was not solved. 

Let us show the Lyapunov stability of Thomson's configurations on a sphere and 
plane relatively to the perturbations of mutual distances between vortices (further 
we shall speak about "relative" stability) for N — 2, . . . , 6. 

The case of N = 2 

The most general kind of motion of two vortices on a sphere (plane) is the uni- 
form rotation around a motionless axis, which passes through the center of vorticity 
and the center of sphere (perpendicular plane) with preservation of the mutual dis- 
tance E] • Under a small perturbation in this system the distance between vortices 
changes a little, but it is also preserved during the motion. Thus, Thomson's con- 
figurations of two vortices are "relatively" Lyapunov stable both on a sphere and on 
a plane. 
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The case of N = 3 



Let us consider a geometric interpretation of motion of three vortices on a sphere, 
that was used in |3J |B] . In Fig. 1 we represent the surface of integral 

D = aiMx + a 2 M 2 + a 3 M 3 

in the space of Mi, M 2 , M 3 , where cij = J- are inverse intensities of vortices, 
and Mi are squared distances between vortices. The triangular area presented in 
Fig. 1 is selected by the conditions Mi > 0. 

We have painted the areas of nonphysical states, for which the triangle inequality 
is not fulfilled over, in black. The curves represent levels of the energy integral. 
Thomson's configuration on this figure corresponds to the point A, where Mi = 
M 2 = M 3 . As we can see from the figure, under small perturbation of Thomson's 
configuration the system starts moving along the level of the energy integral and D, 
close to the point A, which is bounded. Thus, after perturbation the system moves 
in a small neighbourhood of the point A, without leaving it. Analogous arguments 
can be carried out for the case of a plane. Hence, Thomson's configurations of three 
vortices are "relatively" Lyapunov stable both on a sphere and on a plane. 



The case of N = 4, 5, 6; O G [0, 0*) 

To show the "relative" Lyapunov stability of Thomson's configurations for iV = 
4, 5, 6 and for 9 G [0, 9q) let us diagonalize the quadratic part of the Hamilto- 
nian H 2 . Such a diagonalization, as it has been already said above, is carried out 
with the help of Fourier transformation given by the matrix 

T kl = —e~^ kl . 



After this diagonalization and elimination of a degree of freedom, corresponding to 
the motion of vortices as a solid-state configuration, we obtain: 



= ^ cot 9M + q\ + P j + ql) + ^ V3C X7o 1 {PI + ql) ' 

= - cot 6M + q \ +pl + ql) + r W3(2cos^ -l) , + 2 + p2 2 
ii 2n sin 6/q 

rr(6) ,,,2,2,2, 2^ , ^ V5 COS 2 6 - 4 „ 2 
H 2 =~T O (Pi + Q.X+P2 + W + "j — Q (P5 + Qs) + 

An 47r sm o 

r 2 v/2(5cos 2 o -3) ^ 

+ Y, ^ + %+P4 + g 4 ). 
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For 9 G (0, 0q) the quadratic Hamiltonians (JTTJ) are positively defined and, 
consequently, can be used as Lyapunov functions. Thus, Thomson's configurations 
of four, five and six vortices are "relatively" Lyapunov stable on a sphere in areas 
of their linear stability. Proceeding to limit — ^ 0, it is also possible to show a 
"relative" Lyapunov stability on a plane. 



3 The case AT = 7 on a plane 

Thomson's configuration are unstable in linear approximation on a sphere in case 
N = 7. On a plane the linear instability vanishes, but appears in the system of two 
zero frequencies with two-dimensional Jordan cells (the neutral case of a indifferent 
equilibrium). After diagonalization the quadratic Hamiltonian of the system takes 
the form 

h 2 = \{ q l + 4) + y (^ + ll + Pl + qI) + f (P§ + <i + Pl + qI) , (12) 

where a = -, U\ = j^, u 2 = t~, and the values of T and R are assumed to be 
1. It can be seen that this system has two resonances of type 1:1 (variables 3, 4 
and 5, 6) and two zero frequencies with two-dimensional Jordan cells (variables 1 
and 2). The Hamiltonian (|12|) is not positively defined, and, therefore, to study 
its stability it is necessary to investigate expansions of the Hamiltonian of higher 
order. For simplification of the Hamiltonian form we shall carry out the procedure 
of Deprit-Hori normalization (a detailed exposition of the procedure of Deprit-Hori 
normalization with the presence of zero frequencies with two-dimensional Jordan 
cells is presented in Appendix). After normalization the Hamiltonian of the third 
degree will have the form: 

#3 = 6(5 (Pi - Pi)(pI + <&-pI- ll) + (Pi + P2)(P3P4 + qm)) , ( 13 ) 

where b = 0.0672552. H 3 does not depend on the variables p$, q$, pe, q^. Therefore, 
degrees of freedom, corresponding to these variables in approximation of the third 
order, are separated. Further we can consider a system with four degrees of freedom 
only. 

To prove the instability of this system let us use a method of constructions of 
quasihomogeneous truncations of expansion of a Hamiltonian [7 . Thus for a system 
obtained from the initial one by quasihomogeneous truncation, we try to find sets 
of partial solutions as twisted rays ~ rg-sin(o;lnt + a). We have considered basic 
quasihomogeneous truncations of a system of the third degree. However, non of 
them has a partial solution in form of twisted rays. 

The book [7j presents constructions of increasing solutions for systems with 
two degrees of freedoms with resonances 1:1 and two zero frequencies with two- 
dimensional Jordan cells. In the case of resonance 1:1, to construct the increasing 
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solutions we must have the terms of no less than the 4 th degree in variables, corre- 
sponding to the resonance; and in case of two zero frequencies with two-dimensional 
Jordan cells the terms must be no less than of the third degree. Due to particularity 
of the system and a large number of degrees of freedom, there is no such terms in this 
case, H 3 includes only "mixed" terms (depending both on variables, corresponding 
to the resonance 1:1, and on variables, corresponding to zero frequencies). Probably, 
this is the reason why it is impossible to find increasing partial solutions. 

Motion equations of the system of the third degree have an invariant manifold 

P2 = q2 = P3 = 93 = P4 = <?4 = P5 = Q5 = Pe = Qe = . (14) 
The equations of motion on this manifold take the form: 

pi = -aqi , 
qi =0. 

A solution of the given system is a set of particular solutions, which increase linearly 
with time. However, a linearly increasing solution of the truncated system can not 
be extended to an asymptotically increasing solution of the complete system (which 
is possible for solutions in the form of twisted rays). Therefore, we shall consider 
the further expansion of the Hamiltonian up to terms of the fourth degree. 
The normal form of the 4-th degree can be written as 

H 4 = 0.943361gij0ig 2 p 2 + 0.728898W/i 2j R 5 6 - 0.518170W34#56+ 
+ 0.010176i? 2 6 - 0.518170W34W56 + 0.378841^34^56+ 
+ 0.003197^ 2 + 0.513722 W12W34 - 0.120745 J R 12 # 34 + 
+ 0.005339W 3 2 4 - 0.76260^4^34 - 0.559242J¥ 12j R 34 - (15) 

- 0.080602J¥ 34j Ri 2 + 0.353421J¥ 56j R34 + 0.237601Wf 2 - 

- 0.013751iy 56J R 56 - 0.006875W 5 2 6 + 0.769745^i 2 W 56 - 

- 0.030918i4 , 

where 

W12 = P\q2 - P2<?i , Ri2=pl+pl, 

W u = p 3 q 4 - p 4 q 3 , Ru = \ + ll + Pi + » 

W m = p 5 q 6 - p 6 q 5 , R ^ = \ (#3 + 95 + Pi + 9e) ■ 

In the case of H 4 , an analogy with H 3 , we did not manage to find a quasihomogeneous 
truncation of the system which would admit increasing partial solutions in the form 
of twisted rays. 
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Taking into account terms of the fourth degree in the Hamiltonian, equations on 
the manifold (|14|) have the form: 




(16) 

where k = 0.127892. The equations ()16j) have an integral 

k 4 + ^ = const. 
4 2 

Levels of the integral on a phase plane pi, q\ represent closed oval curves of the 
fourth degree with their center in the origin. Hence, motions of the system on the 
manifold ()14j) will be limited. Thus, the addition of terms of the fourth degree in the 
Hamiltonian has resulted in smoothing of linear instability, appearing in analysis of 
the Hamiltonian of the third degree. 

The further normalization is connected with large computing complexity and, 
strictly speaking, the problem of stability for N = 7 remains open. This problem, 
in a sense, is a limiting problem for systems, containing a parameter, which is the 
curvature in the problem under consideration. For the zero value of the parameter a 
linear analysis is insufficient, and we have to take into account the nonlinear terms. 
In favour of instability we can specify the fact that, with adding of small curvature, 
the solutions become unstable in linear approximation. However, this reason cannot 
be considered strict proof of instability. 



4 Conclusion 

Thus, the following generalization of Thomson's theorem on stability of regular N- 
gon configurations of vortices on a sphere is true: 

Theorem 1. Let us consider regular N-gon configurations of vortices on a 
sphere, located on the same latitude 9q\ 

1. For N > 7, they are unstable with respect to deformations; 

2. For N = 6, 5 and 4, they are Lyapunov stable with respect to deformations 
for 9 e (0, 0q(N)), and unstable for 9 e (9q(N), |). The limiting latitudes of the 
stability 9q (N) are determined by formulas 

3. For N = 2, 3, they are Lyapunov stable with respect to deformations for all 9®. 
Remark 1. The problem of stability at 9q = 9q at N = 4, 5 and 6 requires 

separate analysis with respect to nonlinear terms and remains unsolved to the present 
day. 
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As for a plane, we have the following: 

Theorem 2. Regular N-gon configurations of vortices on a plane are Lyapunov 
stable with respect to deformations for N < 7, and unstable for N > 8. The problem 
of stability for N = 7 remains unsolved. 



5 Appendix 

In this work we use the method of Deprit-Hori normalization for construction of 
normal forms of Hamiltonian near a fixed point. A detailed exposition of this method 
in case of absence of Jordan cells can be found, for example, in Here we shall 
generalize this algorithm to the case of the presence of two-dimensional Jordan cells. 
To determine a normal form we shall apply the transformation close to identical, 
defined by a Hamiltonian W. The Hamilton form of transformation ensures its 
canonical properties. The initial Hamiltonian H, a normalized Hamiltonian K and 
the function W are expanded into a series in terms of a small parameter near a fixed 
point 

oo oo oo 

n=0 n=0 n=l 

Here all functions H n , K n and W n are polynomials of n degree with respect to 
phase variables. A normalized Hamiltonian up to k degree is obtained by successive 
solution of operator equations 



K = H , 

n-2 



K n = H n + Y] iL^xH^^x + r7#i n ~ m) 

i L In — m)l 



m=0 



(17) 

+ L n H 0l n > 1 , 



where 



3=0 

for n = 0, 1, . . . , k. 

Thus, on the n-th step, the normalization procedure is reduced to the solution 
of the equation 

K n = H* n + L n H , (18) 

where H^ n is a known function, depending on functions already found on the pre- 
vious steps of normalization. 

As we consider the normalization near a fixed point, it is obvious that the ex- 
pansion of H will begin with the second order terms. Hence, the corresponding 
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expansions of K and W will begin with the second and the third order terms. The 
equation (|18j) can be rewritten in the form 

K n = H k n n + VW ni (19) 

where Vf = {H 2 , /}. 

While solving the equation (|T9*|) . the function W n is selected so that all terms in 
the right part of the equation (fT9"j) are excluded, except those, which could not be 
presented as M = T>f, where / is some function. K n will consist of these remained 
terms. 



5.1 Absence of Jordan cells 

Let a quadratic part of the Hamiltonian H 2 lack Jordan cells. Then after diagonal- 
ization and complex change 

1 t 

z = p + iq, z = p-iq, {z k ,Zi\ = —d k i 
the quadratic part of the Hamiltonian H 2 has the form: 

% 

The action of an operator T> on some monomial Mi m = Ai jTn z l z m is defined by 
expression 

VM l;m = -u>{m- l)M l>m = C ltm M l>m , (20) 

where u> is a vector of frequencies. The operator T> does not change the degree of 
monomial, but only multiplies it by some constant. 

Thus, in the right part of (|17|) we are unable to eliminate only those monomials, 
for which the following equality 

Ci >m =-^{m-l) = (21) 

is fulfilled. A normalizing function and a normalized Hamiltonian will take the form: 

M, m 



K n = V M iirn , W n = V r< 



l,m l.m 



l.m 



where Mi m are monomials included in H kn . 
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5.2 A case of two-dimensional Jordan cells 

Let H2 have k zero frequencies with two-dimensional Jordan cells, then H2 can be 
presented as 

k N 



#2 = XX+ £ ^i+pI 



i=l i=k+l 



After a complexification with respect to variables qk+i, • • • , QV, Pk+i, ■ ■ ■ , Pn we 
obtain 



k N-k k 

i=l i=l i=l 

In this case, the operator D has the form 

D = D' + ^ 2gi — , (22) 

where the operator D' acts on a space of variables zi, . . . , %_ fe , zi, . . . , %-fc> and 
has the same properties, as the operator D in Section 5.1 (i.e. it does not change 
the degrees of monomials). In this case, the operator D does not already preserve 
the form of monomials. Let us consider its action on some monomial 

where 



z l ■ ■ ■ z N-k 2 1 " • • • Z N-k 



11^7 ^Sl ^"N-k — ri — ""JV— fe 



Applying the operator 2) to Mj jTTl we obtain 

k 



®M hm = C s , r M l<m + Mi m J22U-, C s<r = &M' a r = ^uj(r-s). (23) 

^ Pi 2% 



i=i 
Let 



l,m, l.m l,m 

then the equation (|19|) for coefficients of monomial M( ra will take the form 

= A tl + B i% C *>r + £ Btl, m - ei 2(h + 1) • (24) 

i=l 

Here ei is unit vectors of an integer lattice. From (|24|) it is clear that it is impossible 
to set Kj^t e Q ua l to zero, if both conditions are fulfilled simultaneously: 
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1) C s r = i. e. M' s r is a resonance monomial; 

2) rrii = 0, % = 1, ...,k, thus the sum in (|24j) is absent, since for negative 
components l,m do not exist. 

The monomials which satisfy these two conditions are resonance monomials. 
However, as we shall see further, there are also other resonance monomials. 

Let us obtain now a normalizing function, which will help us to exclude nonres- 
onance monomials. Let us consider separately 2 cases: 

I. M' „ is not resonance: 

II. M' sr is resonance, but rrii ^ 0. 

Case I. Let us consider an integer lattice l\, . . . , 1^. For each monomial Mi )Tn and 
function wi tTn , which excludes the latter, we shall put in the correspondence a point 
of this lattice. Let us denote L-plane as the hyperplane cutting off a fc-dimensional 
angle in Z-space, determined by the equation 



8=1 



Let H^ n contain a monomial Mj m . The given monomial can be excluded from the 



right part of (fTT?|) with the help of function Wi tTn = — C ,' ,T " . Under this elimination, 
as we can see from (|2*Hj) . in the right part of (jl9j) we get the sum of monomials 

'M £ _ ei>m+ei . (25) 



k 



2k 



=1 



Thus, we pass from the point I to the neighbouring points I — e^, i = 1, . . . , k. Such 
passage is always made from L-plane into (L — l)-plane, and the monomials are 

multiplied by I —-M^- I . Using the same method we exclude originating monomials 



Cs, 

and move along a lattice in the direction to the origin. In the general case an 
arbitrary point of the lattice can be reached by several paths, and at the same time 
appearing monomials are summarized. When any coordinate becomes equal to zero 
at such a motion, the further motion takes place in remaining variables only. In 
Fig. 2 the lattice and such passages at k — 2 are presented. A sequence of the 
passages constructed in such a way breaks on a point 1 = 0. 

Adding normalizing functions wi iTn by all points of a rectangle, we obtain a 
function Wi )7ni that excludes a monomial Mi m : 

w ^ = -c£ E M '--"m^ij- v czJ f '* (28) 

where m' = m+(l — l'), fin is a number of paths from point I into point I' moving 
in direction of origin. This number can be calculated with the help of recurrence 
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relation 

k 

fi',i = y] fi'+et,ip- - <%zj ■ 

i=l 

Besides the successive elimination of monomials of the first type with the help of 
functions (|26j) . they can also be eliminated with the help of functions u>z jTn ,. For 
this method it is neccessary to take such order of enumeration of Ui m that with 
elimination of the next monomial in the right part of (|19|). the excluded earlier 
monomials do not appeared again. Since the elimination of monomial Mi m with 
the help of wi )Tn we pass from L-plane into (L — l)-plane, such order is a sequential 
enumeration of all monomials with maximum L = n, then with L smaller on one 
unit, etc. 

Case II. In this case an action of the operator T> over monomial Mi y7n has the 
form 

k 

VM hm = 2Z i M l _ ei)Tn+e4 . (27) 
i=i 

Thus, it is possible to exclude the monomial Mi >m with the help of the functions 

< = -2ftTT) Mi+ — ■ (28) 

And in the right part of (J19)) the sum 

I, 

(-- 



h + 1 



M l+ei - e.j ,m — ei-\-ej 

(29) 



is obtained. As we can see from the interpretation of an integer lattice, with such 
elimination we pass from the point / in points I + e$ — ej, j = 1, . . . , k, j ^ i. 
Thus there occures a motion in one L-plane. 

Successively applying transformation with the functions (|2*%|) for various i we ob- 
tain some chain of monomials. Selecting the coefficients for excluding functions (|2*H|) 
we can reduce the whole chain to any of its monomials. Thus, if we can exclude 
only one of monomials from a chain without emerging of new terms, the whole 
chain can be excluded from the right part of ()19|). The monomial can be completely 
excluded, when all terms in the sum (j2l?|l equal zero, i. e. when for some of % we 
have lj = 0, j ^ i and m, > 0. We can reduce an arbitrary monomial to such form 
(and exclude it completely), if there exists at least one i, such that 

mi>Y^ lj . (30) 
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We can prove it with the use of an invariance of the inequality (|3U|) with respect to 
the transformations (J28j) . If for the monomial Mt, m the inequalities 

rrii < lj , % — 1, . . . , k (31) 

are fulfilled, the given monomial cannot be excluded completely from the right part. 
We can only transfer it into any of monomials from the chain. Thus monomials, 
satisfying (}3"T]) are resonance monomials. 

Let us consider the case, when for some i the inequality ()30|) is fulfilled, i. e. the 
monomial can be excluded completely. Let us construct a function Wi^ m while h 
will exclude the monomial. After an elimination of this monomial with the help of 
the functions (J28j) we pass from point I into points I + — e^, j ' ^ i. On L-plane 
the motion takes place in direction to an axis U. At k = 3 such a plane and the 
sequence of passages are represented in Fig. 3. For each passage of monomial from 
some point into another one, the monomial is multiplied by — As in case I, 
summing up wi jTn by (A; — l)-dimensional rectangle, we obtain 

wu = E- E E ■■■E(- 1 >' H ' +1 np "f. 2 r/'+ 1 ) Mf+ "-'-- (32) 

where 

^ = h+^ilj-l'j), m'j = nij + ilj-lj), i = l, ...,k f v ,i = /i/_ e . +e ., i (l-^ iij ) . 

Thus, in the case of k zero frequencies with two-dimensional Jordan cells normal 
form of the Hamiltonian include the monomials Mj ;Tri = M' s r p 1 ^ . . -P^q™ 1 ■ ■ ■ q™ k 
such that M' s ^ r are resonance and the following inequality 

m i < lj , for V i = 1 , . . . , k 

is fulfilled. 
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